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New upper bounds on direct detection rates have recently been presented by a number of experi- 
mental collaborations working on searches for WIMPs. In this paper we analyze how the constraints 
on relic neutralinos which can be derived from these results is affected by the uncertainties in the 
distribution function of WIMPs in the halo. Various different categories of velocity distribution 
functions are considered, and the ensuing implications for supersymmetric configurations derived. 
We conservatively conclude that current experimental data do not constrain neutralinos of small 
mass (below 50 GeV). 

PACS numbers: 95.35.+d,98.35.Gi,11.30.Pb,12.60. Jv,95.30.Cq 



I. INTRODUCTION 

In Rcfs. Q, 1^, Q we have discussed the cosmologi- 
cal properties of light neutraUnos (i.e. neutraUnos with 
a mass in the range 6 GeV < < 50 GeV), which 
originate in supersymmetric schemes where gaugino-mass 
unification is not assumed. Actually, the most remark- 
able features occur for neutralinos in the mass range: 7 



for relic neutralinos), one recovers the expression: 



dR 
dE^ 
where: 



I[Vmin) = I d-'w 

'l«>'l'mi„ ^ 



(2) 



In the previous formulae, notations are: Nt is the num- 
GeV < mx < 25 GeV. Namely, for relic neutralinos with ^er of target nuclei per unit mass, is the WIMP 



these masses, direct and indirect detection rates are con- 
siderably high, and at the level of present experimental 
sensitivities. Furthermore, the range of the predicted val- 
ues for the rates is quite narrow, at variance with what 
happens for neutralinos of higher masses, where the ex- 
pected rates are spread over decades. 

The properties of light neutralinos with respect to 
WIMP direct measurements were analyzed in Refs. |^|^. 
Since then, some new results and/or analyses of previ- 
ous data from experiments of WIMP direct searches have 
appeared [BiSSB- In the present paper, we examine 
whether these new data put some constraints on the relic 
neutralinos of light masses. 

Let us recall that the differential event rate dR/dEn 
{Er being the nuclear recoil energy) measured in WIMP 
direct searches is a convolution of the WIMP-nuclcus 
cross section with the WIMP phase-space distribution 
function of WIMPs, evaluated at the Earth location. By 
assuming that: i) in this phase-space distribution func- 
tion, spatial and velocity dependence factorize, and 2) 
coherent interactions dominate over incoherent ones in 
the WIMP-nucleus scattering (which is usually the case 



mass, TTiTv is the nucleus mass, /ii is the YflMP-nucleon 
reduced mass, A the nuclear mass number, crg""[^'™"^ 
is the WIMP-nucleon coherent cross section, F{Eii) is 
the nuclear form factor, ^ is the fraction of the mass 
density of the WIMP in terms of the total local den- 
sity for non-baryonic dark matter po (i-e-: S, ~ Pw/po)i 
fEs{w) and w denote the velocity distribution function 
(DF) and WIMP velocity in the Earth frame, respec- 
tively {w = \w\). It is natural to define the velocity dis- 
tribution function in the Galactic rest frame f{v), where 
V = w -\- V(Q, Vi^ being the Earth velocity in the Galactic 
rest frame. The Earth frame velocity DF is then obtained 
by means of the transformation: /es(w') = f{w -\- '?©). 
It is implicitly understood that the velocity DF f{v) is 
truncated at a maximal escape velocity Wcsc, since the 
gravitational field of the Galaxy cannot bound arbitrar- 
ily fast WIMPs. The value we adopt here is: Vcsc — 650 
km sec~^ |^, although we will comment on the effect 
of a lower value, which we will set at Vcsc = 450 km 
sec~^ Q. Finally, the quantity v^^in appearing in Eq. (0) 
defines the minimal Earth-frame WIMP velocity which 
contributes to a given recoil energy Ej^: 



fmin = [mNER/{2n\)]^/\ 



(3) 



A: Spherical puM., isotropic velocity dispersion 

AO Isothermal sphere 

Al Evans' logarithmic 

A2 Evans' power-law 

A5 NFW [16] 



B: Spherical pdm, non— isotropic velocity dispersion 

Bl Evans' logarithmic [14j 

B2 Evans' power-law 

B5 NFW [16] 



C: Axisymmetric puM 

C2 Evans' logarithmic 

C3 Evans' power-law 



D: Triaxial ^dm fl7] 

Dl Earth on major axis, radial anisotropy 



TABLE I: Summary of the galactic halo models considered 
in our analysis. The label shown in the first column is used 
throughout the text to indicate each model in a unique way 
and corresponds to the classification introduced in Ref. ^0|. 
For details on the models and proper definitions, see Ref. |lO|. 
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FIG. 1: Function 2r(iimiii) for an isothermal sphere (model 
AO of Table EJ and u^sc = 650 km sec . The solid curves 
(from top to bottom, as seen on the extreme left of the plot) 
refer to the three values of vq (and ensuing po) given in Table 
im vq — 170 km sec~^ (top), vo — 220 km sec~^ (medium), 
vq — 270 km sec~^ (bottom). The dashed line shows the 
modification of the median isothermal case when Dcsc = 450 
km sec~^. 



where fiA is the WIMP-nucleus reduced mass. 

Eqs. lU and © are the basis for deriving information 
on the quantity C'^'scria'™"^ from the measurements on 
the differential rate dR/dEn. However, this procedure 
imphes the use of a specific WIMP distribution function, 
which determines both the value of the local dark matter 
density po and the shape of the velocity DF f{v). 

In the present paper we first discuss how upper bounds 
on '^crg™j?|'™"' , derived from experimental upper limits on 
dR/dEn, depend on the large uncertainties affecting the 
WIMP distribution functions. We then discuss what is 
the relevance of these upper bounds on C'^gcaia™"^ 
light relic neutralinos. 

II. WIMP DISTRIBUTION FUNCTIONS 

In our analysis we consider a subset of the large sam- 
ple of galactic halo models which were studied in detail 
in Ref. 0| . Following Ref. ^3 ' classify the DFs into 
four categories, depending on the symmetry properties 
of the matter density (or the corresponding gravitational 
potential) and of the velocity dependence: A) spheri- 
cally symmetric matter density pdm with isotropic veloc- 
ity dispersion, B) spherically symmetric matter density 
with non-isotropic velocity dis pers ion, C) axisymmetric 
models, D) triaxial models [lllll2j 



For each category, different specific models are identi- 
fied. The models considered in the present analysis are 
listed in Table For a thorough definition of the differ- 
ent models and of the values of their intrinsic parameters, 
and for a detailed description of theoretical technicalities, 
we refer to Ref. jl^l. Here we just remind that for each 
model we calculate, either analytically (when possible) or 
numerically, the velocity DF which accompanies a given 
matter density distribution. For the spherically symmet- 
ric and isotropic models of class A, the velocity DF is 
obtained by solving the Eddington equation [Hill. For 
the spherically symmetric and non-isotropic models of 
class B we assume the anisotropy to be described in terms 
of the Osipkov-Merrit parameter /3 0,0,0] which de- 
fines the degree of anisotropy (we fix /3 = 0.4): in this 
case, the velocity DF can be obtained by a generaliza- 
tion of the Eddington method 0, Axisymmetric 
models of class C allow the presence of a definite angular 
momentum. We choose them as a direct generalization 
of some of the models of class A: for these models, an- 
alytical solutions to the relevant generalized Eddington 
equation may be found. In this case, we also allow for a 
maximal co-rotation or counter-rotation of the galactic 
halo. Finally, class D presents a specific triaxial model. 

Each halo model is constrained by a number of observa- 
tional inputs [lol | : i) properties of the galactic rotational 
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vo = 170 km sec ^ 


vq = 220 km sec ^ 


t;o ~ 270 km sec ^ 


Model 


po (GeVcm-3) 


po (GeVcm-^^) 


po (GeVcm-^) 


AO 


0.18 


0.30 


0.71 


Al , Bl 


0.20 


0.34 


1.07 


A2 , B2 


0.24 


0.41 


1.33 


A5 , B5 


0.20 


0.33 


1.11 


C2 


0.67 


1.11 


1.68 


C3 


0.66 


1.10 


1.66 


Dl 


0.50 


0.84 


1.27 



TABLE II: Values of the dark matter local density po corresponding to the three different values of the local rotational velocity 
Vo and obtained from the constraints on the amount of non-halo component and on the flatness of the galactic rotational curve, 
for the different halo models of Table |I] For the models of class A and B, the values of po are the minimal ones for vo = 170 
km sec~^ and vo = 220 km sec"'^ (i.e. corresponding to a minimal halo contribution), while for vo = 270 km sec~^ the values 
of Po are the maximal ones (referring to a maximal halo). For models of class C and D, the value of po is always the maximal 
one. The axisymmetric models of class C are not affected by the inclusion of a co-rotation or counter-rotation effect. 
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tribution of the halo to the rotational curve is minimized, 
the value of po is either increased or reduced, respectively. 
This is discussed in detail in Ref. ^3 a-^d is manifest also 
in Table im where values of po for each halo model and for 
the three representative values of vq are shown. The dif- 
ference in the values of po depending on the assumption of 
a minimal or a maximal halo contribution is given in Ref. 
[lo|. For instance, for the isothermal sphere po falls in 
the range (0.18,0.28) GeV cm'^ for vn = 170 km sec"!, 
(0.30,0.47) GeV cm-3 for vq ^ 220 km sec^^, and 
(0.45,0.71) GeVcm-3 for = 270 km sec" i. We 
therefore expect a significant variability in the extracted 
upper limits in direct detection experiments also for 
the simple and generally used isothermal sphere model. 
We notice that the standard reference choice of po = 
0.30 GeV cm^"^ for vq ~ 220 km sec~^ refers to the case 
of a minimal halo. 



FIG. 2: Values of Umin as a function of the nuclear recoil 
energy iJ_R, for a Ge detector. The different curves refer to 
WIMP masses of: 10, 20, 50, 100, 200 GeV and 1 TeV, from 
top to bottom. 



curve, namely the range of the allowed values for the local 
rotational velocity, 170 km sec^^ < < 270 km sec^^ 
[iolEH . and the amount of flatness of the rotational curve 
at large distances from the galactic center, and ii) the 
maximal amount of non-halo components in the Galaxy, 
Mvis (i.e. the disk, the bulge, etc.). These constraints 
determine first of all the value of the local dark matter 
density which is a relevant parameter in the direct 
detection rate. Depending on whether one allows for a 
maximal halo [i.e. the contribution of the non-halo com- 
ponents is minimized) or, the other way around, the con- 



The choice we make here for the values of po associated 
to each representative value of wq is the following: for 
Vq = 170 km sec^^ (which correspond to its 95% C.L. 
lower bound) we adopt the case of a minimal halo, in 
order to determine the set of less constraining upper- 
limits on ^CTscaia'r ' clcar from Eq. For vq = 
270 km sec~^ (which correspond to its 95% C.L. upper 
bound) we instead adopt a maximal halo: in this case 
we will determine the most constraining upper-limits on 
^''^scaiar™"'' ■ Central (and reference) 
value Vq = 220 km sec""'^, we adopt a minimal halo, which 
reproduces the standard choice po = 0.30 GeV cm^'^ for 
the isothermal sphere. These considerations apply to all 
models of class A and B. For the models of class C and 
D, for which we can rely only on analytical solutions for 
the velocity DF, we are forced to use always the case 
of a maximal halo: in fact, analytical solutions of class 
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FIG. 3: Function X(uniin) for all the galactic model of Table HJ other than the isothermal sphere, for v^sc = 650 km sec~^. 
The label which identifies the model is written in the bottom-left corner of each panel. Notations are as in Fig. Q In panels 
(g) and (h), which correspond to axisymmetric models, the dotted and dashed lines refer to maximal galactic co-rotation and 
counter-rotation, respectively. In panel (h), the long-dashed line shows the modification of the vo = 270 km sec~^ case when 
fcsc = 450 km sec~^. 



C and D are actually obtained only for a maximal halo 
contribution [lfll |. 

Let us turn now to the direct effect of the velocity DF 
on the detection rate, which is studied here in terms of 
the relevant hmction I{v,^{^) of Eq. ISJ. Fig. shows 
I{v-min) for thc isothcrmal halo and for the three values 
of vq listed in Table |nj We sec that for low values of Vmin 
the larger contribution to the detection rate occurs when 
Vq is smaller, since for smaller rotational velocities the 
velocity dispersion of the isothermal Maxwellian distri- 



bution is also smaller, and in turn this enhances the av- 
erage inverse velocity, which is related to the definition of 
I{vmin)'- {^/w) = 2"(0). This can be analytically under- 
stood by remembering that for a pure isothcrmal sphere 
and a maximal halo the velocity distribution function 
is just an isotropic Maxwellian with velocity dispersion 
given by wq: 

/^(«) = [7r^;2]-3/2cxp(-^;V«o'); (4) 
in this case the function reads (in the limit v^sc — * 
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oo) HI: 

^(wmin) = T7^[erf(a;niin + ri) - erf(x„iin - rj)] , (5) 

where Xmin = I'min/i^o ^nd 1] = Vi^/vq. From Eq.© we 
see that for small values of Umin the larger J{vniin) occurs 
for smaller vq because of the inverse law dependence. 

On the contrary, for large values of Wmin, the almost- 
exponential tail in I(t'min) is more severe when vq is 
small, and therefore the behaviour of X(i;min) with re- 
spect to vq is the opposite. This again is understood 
from the simple expression of Eq.©. The regime we are 
considering (wmin ^ I'e) asymptotically can be studied as 
the limit 77 — > in Eq. jSJl, which gives: 

I{vmin) = 2[7:vl]-^/^ cxpi-vl^Jvl) . (6) 

This shows the discussed behaviour as a function of vq. 
The tail, due to the presence of a non vanishing 77 in Eq. 
Q , is less severe than the one in Eq. ^ but nevertheless 
it follows the same behaviour. 

Also the value of the escape velocity is relevant in the 
large Wmin tail of the function X{vmin)- The results pre- 
sented so far in Fig. ^ are obtained for a value of the es- 
cape velocity (in the galactic frame) v^sc = 650 km sec~^. 
However, values as low as Vesc = 450 km sec~^ have also 
been considered . A lower escape velocity implies a cut 
in the high Wmin tail of T('ymin) |22| . This effect is shown 
in Fig. ^for the central case vo = 220 km sec^^. 

The discussion on the behaviour of X(wniin) has direct 
impact on the direct detection rate, since Wmin is directly 
related to the recoil energy. Eq. © implies that very 
light WIMPs can produce recoil energies in the tens of 
keV range only if they possess large velocities. In this 
case, the detection rate for such light WIMPs will be 
mostly determined by the almost-exponential tail in the 
function T(wmin) discussed above. On the contrary, heavy 
WIMPs can produce recoil in the same tens of keV range 
by possessing much lower velocities: they will be there- 
fore more sensitive also to the low Wmin part of the func- 
tion X(v„ii,i). The quantitative connection between Wmin 
and Eji for a Ge nucleus and different WIMP masses is 
given in Fig. |2| 

Finally, Fig. |21 shows the function Z(wi„in) for all the 
other halo models listed in TableQ] As for the symmetric 
and isotropic models, we see that in the case of a power- 
law behaviour of the gravitational potential (model A2) 
or for the NEW density profile (model A5) the large Wmin 
tail is less suppressed, mainly for low vq: this comes along 
with a larger predicted detection rate and it will translate 
into a more constraining upper limit for WIMPs lighter 
than a few tens of GeV. In the case of anisotropic models, 
we notice that the most direct anisotropic generalization 
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FIG. 4: The solid lines show the upper limit on the quantity 
Cc^scTiar"' as a function of the WIMP mass for the CDMS 
detector and for an isothermal sphere (model AO of Table 0. 
The curves refer to the three values of vq (and corresponding 
po) given in Table UTI vq — 170 km sec~^ (top), vq = 220 km 
sec~^ (medium), vo — 270 km sec~^ (bottom). The colored 
region shows the values of ef^(""'^''=°") fQj- neutralino dark mat- 

o ^ scalar 

ter, obtained in a scan of the minimal supersymmetric model 
defined in Refs. The funnel for low neutralino 

masses (below 50 GeV) corresponds to supersymmetric mod- 
els without gaugino-mass unification. 

of the isothermal sphere, which is a cored spherical distri- 
bution with anisotropic velocity dispersion (model Bl) is 
the one which is more suppressed at large Wmin: this has 
the effect of reducing the sensitivity of the detector to 
light WIMPs. On the contrary, the axisymmetric model 
with a power-law gravitational potential C3 is the one 
with the highest tail in the function X(v^in), and there- 
fore more sensitive in constraining light WIMPs. We also 
notice that for this type of models, which possess an cn- 
haced Wmin tail, the effect of a lower escape velocity is 
more dramatic: panel (h) of Fig. O shows the sizeable 
reduction of X(?;min) at large Umin for an escape velocity 
Vcsc = 450 km sec~^ in the case of model C3 and vq ~ 270 
km sec~^. 

The local matter density values po of Table |n] and the 
results of Fig. ^ and Fig. |2| are the key elements which 
will be used in the next section to determine upper limits 
on the WIMP-nucleon scattering cross sections. 

III. RESULTS AND CONCLUSIONS 

In Refs. 0, H, H 13 upper limits on ^cTsc^iIi"'^^ ^^'^ o^)- 
taincd by using a standard isothermal distribution tunc- 
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FIG. 5: Ratio of the upper limits of Fig. |1] obtained for an 
isothermal sphere (model AO). The upper curve is the ratio 
between the upper and the central curves in Fig. |31 The lower 
curve is the ratio between the lower and central curves 



tion with density parameter po — 0.3 GeV cm~'^ and 
vq = 220 km sec~^. Here, we analyze this class of exper- 
imental data employing the sample of distribution func- 
tions discussed in the previous section. Our analysis and 
discussion is performed in terms of the experimental data 
of Ref. (CDMS), since these turn out to be the most 
constraining ones. For very light WIMPs, with masses 
below 10 GcV, also the rcsuhs of Ref. (ZEPLIN) may 
play a role, and we will add them in our final discussion. 

In order to obtain the exclusion-plot from the CDMS 
data, we extract the neutralino-nucleon cross section 
that yields a number of events compatible with zero be- 
tween 10 keV and 64 keV with an effective exposure 
(MT)off = 19.4 kg day, which corresponds to the effec- 
tive exposure for = 60 GeV quoted in Ref. |^. By 
considering a Poissonian fluctuation of the expected rate, 
we assume the upper bound of 2.3 events at 90% C.L. In 
the calculation of the expected rate, we use Helm nuclear 
form factors and a bolometric quenching factor equal to 
1, as quoted by the experimental Collaboration [8j. Our 
procedure for extracting the exclusion-plot is less refined 
than the one adopted in Ref. 0|, since it neglects the 
dependence of (MT)efi on the WIMP mass and does not 
take into account the part of the CDMS spectrum be- 
tween 64 and 100 keV (which can be included follow- 
ing for instance the statistical procedure of Ref. [13 )■ 
Nevertheless our procedure allows us to reproduce to a 
good degree of precision the CDMS limit for the standard 
isothermal distribution, as shown by the central curve in 
Fig. 21 when this is compared with the upper bound dis- 



played in Fig. 39 of Ref. [3. Therefore, in our analysis 
we adopt our simpler procedure; we have checked that 
adding a proper treatment of the efficiency and adopting 
the statistical procedure of Ref. yields quite similar 
results. 

The upper limits for the isothermal sphere (model AO) 
are shown in Fig. 0] for the three representative values of 
Vq and the corresponding choices for the local dark matter 
density po, as quoted in Table ^1 As already mentioned, 
the central curve corresponds to the reference case of uo = 
220 km sec~^ with po = 0.3 GeV cm^'^. The upper and 
lower curves are instead obtained for vq = 170 km sec~^ 
with Po = 0.18 GcV cm~'^ and vq = 270 km sec"""^ with 
Po ~ 0.71 GcV cm~'^, respectively. An important effect 
is obviously due to the different values of po which are 
associated to the different values of vq, as discussed in 
the previous Section and in Ref. However, the dif- 

ference in the function Xivmin) is quite relevant in the de- 
termination of the upper limits, especially at low WIMP 
masses. In order to appreciate the difference in the ex- 
clusion plots, we show in Fig. |Slthc ratios of the upper 
limits obtained with vq = 170 km sec~^ (lower curve) 
and Vq ~ 270 km sec~^ (upper curve) with respect to 
the central vq = 220 km sec^^ case. The dashed hori- 
zontal lines show the ratios of the corresponding values 
of pq. Wc can notice that at low WIMP masses the dif- 
ference in the exclusion plots is very large, much larger 
than the naive ratio of the corresponding po's. This is a 
consequence of the sizable difference in I(urnin) for large 
Vmin; which is the regime relevant for light WIMPs, as 
discussed before. The steep behaviour of the ratios in 
Fig. Ois a consequence of the fact that the sensitivity of 
direct detection to very low WIMP masses (below about 
10 GcV) rapidly vanishes. 

On the contrary, at large WIMP masses, the differ- 
ence in the exclusion plots is much close to what one 
would expect on the basis of the difference in the po val- 
ues. This is clear from our previous analysis, since for 
large WIMP masses the relevant range of Vmin is in the 
100-300 km sec~^ range, which is where the difference 
in Z(f,nin) is small. We can also notice that, for large 
WIMP masses, it is even possible to revert the value of 
the ratio of the exclusion plots naively obtained by the 
ratio of the different po's: this is a consequence of the be- 
haviour of X(uinin) at small Umin discusscd in the previous 
Section. 

Fig. 2] represents the maximal variability which oc- 
curs for the isothermal sphere: this quantifies the as- 
trophysical uncertainty connected to this halo model. 
Confronting the upper limits with the results obtained 
for light neutralinos in supersymmetric models with- 
out gaugino-mass universality 0, Q, IE Bi we can 
see that while all the configuration in the mass range 
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FIG. 6: The solid lines show the upper limit on the quantity ^''"^"^11'™"^ 



as a function of the WIMP mass for the CDMS 



detector and for the different galactic models of Table HJ other than the isothermal sphere. The label which identifies the model 
is written in the bottom-left corner of each panel. Notations are as in Fig. |1] In panels (g) and (h), which correspond to 
axisymmetric models, the dotted and dashed lines refer to maximal galactic co-rotation and counter-rotation, respectively. 



15 (8) GcV < < 25 GcV arc excluded for the central 
(upper) values of vq , only a small fraction arc eliminated 
when Wo assumes its lovirer bound value. Therefore, the 
conservative attitude which has to be taken when set- 
ting limits makes us to conclude that for the isother- 
mal sphere, direct detection only mildly constrains the 
light neutralino sector of supersymmetric models without 
gaugino-mass universality 0, 13, H, ^ , in the 20-40 GeV 
mass range. Clearly the variation on the upper limits due 
to the difference in the halo properties has consequences 
also on the exploration of the supersymemtric parame- 



ter space for heavier neutralinos, as is shown in Fig. 01 
This is relevant also to gaugino-mass universal models, 
for which the lower bound on the neutralino mass exceeds 
50 GeV 0. 

The results for the other galactic halo models is shown 
in Fig. El The differences in the upper limits can be un- 
derstood on the basis of the discussion on the isothermal 
sphere and on the properties of for the different 

models presented in the previous Section. We notice that 
some models, like C3, Dl and B5 are more constraining, 
while in the case of models like Al and Bl the limits 



8 



A, Bottino, F, Donato, N, Fornengo, S, Scopel {2005) 




5 10 50 100 500 

n\ (GeV) 

FIG. 7: The solid lines show the summary of our analysis on 
the upper limit on the quantity ^''"scaia™"^ ^ function of the 
WIMP mass for the CDMS detector and for Vbbc = 650 
km sec^^. The median line refers to the standard isother- 
mal sphere with vo ~ 220 km sec^^ and po ~ 0.3 GeVcm"^ 
(model AO). The upper and lower curves show the two ex- 
tremes obtained in the analysis and refer to model Bl with 
Vo = 170 km sec~^ (upper solid line) and model C3 with 
Vo — 270 km sec~^ (lower solid line). The dashed line refers 
to model C3 with maximal counter-rotation of the galactic 
halo. The dotted lines show the ZEPLIN I limits obtained for 
the same galactic models. The long-dashed lines show the 
upper limits for CDMS in the case of a lower escape velocity 
Vesc = 450 km sec~^: the upper line refers to model AO, the 
lower one to model C3. For model Bl, the limit coincides 
with the corresponding solid line. The colored region shows 
the values of t(j(""'^'<=°") fQj- neutralino dark matter, obtained 
in a scan of the minimal supersymmetric model defined in 
Refs. Q, QIHB' The funnel for low neutralino masses (be- 
low 50 GeV) corresponds to supersymmetric models without 
gaugino-mass unification. 



imposed by direct detection are relatively less severe. 

Finally, we report in Fig. |7|tlic summary of our analy- 
sis: together with the standard central isothermal sphere, 
we show the more (C3) and less (Bl) constraining models 
we obtain. From the analysis of this figure, we conser- 
vatively conclude that from direct detection experiments 
there is currently no constraint on the light neutralino 
sector of supersymmetric models without gaugino univer- 
sality 0, 0, y, Q ■ Should the local value of the rotational 
velocity be on its high range (close to dq = 270 km sec^^) 
direct detection could be able to set stringent limits on 
these supersymmetric configuration: all the mass range 
above 7-8 GeV (depending on the actual halo model) 
and below 25 GeV would be excluded. Notice that. 



would this be the case, also the local density po would be 
large (above 0.7 GeV cm~^, as discussed in Ref. 0| and 
shown in Table HHl : in this case the neutralino configura- 
tions below 7-8 GcV, which arc not constrained by direct 
detection, would be completely excluded by antiproton 
searches However, due to astrophysical uncertain- 
ties which affect the different detection rates, currently 
it is not yet possible to set absolute limits, neither from 
indirect detection techniques 0, or, as shown in the 
present analysis, by direct detection. 

Fig. [7| also shows the effect of a lower escape veloc- 
ity. As discussed in the previous Section, this implies 
a cut in the high Wmin tail of I(t;niin): this turns into a 
weaker sensitivity of direct detection to low-mass neu- 
tralinos. The effect is especially manifest for the most 
stringent models, like model C3 with vq = 270 km sec~^. 
For Vcsc = 450 km sec~^, all neutralino models below 
9 GeV are not constrained even for C3 model. For the 
AO model with vq = 220 km sec~^, there is also a size- 
able difference for light WIMPs, although this is not rel- 
evant for the neutralino configurations. Finally, in the 
case of model Bl with vq = 170 km sec~^, the lower es- 
cape velocity does not produce differences, since in the 
high velocity tail I{v-a_-nn) was already depressed even for 
fosc = 650 km sec~^, as can be seen in Fig. |31 

For completeness, FigO also shows the upper limits 
we obtain, for the isothermal sphere and for the two ex- 
treme cases, for the ZEPLIN detector 0]. The limits we 
obtain for the isothermal sphere arc slightly higher at low 
WIMP masses than the ones quoted by the experimental 
Collaboration Q . We trace this effect to some differences 
in the analysis of the data (we do not make use of the 
"light response matrix" discussed in Ref. Q , since we do 
not have it at our disposal). Fig. |3 shows that for very 
low WIMP masses ZEPLIN could be slightly more sensi- 
tive than CDMS. Nevertheless, even lowering by a factor 
of 2 the upper limits we obtain for ZEPLIN, our conclu- 
sions on the limits imposed to light neutralinos remain 
unchanged. 

Finally, we wish to remind that an annual modula- 
tion effect in direct detection has been observed by the 
DAM A Collaboration over seven years [i^. This result, 
when interpreted in terms of scalar WIMP-nucleus inter- 
actions, leads to an allowed region in the plane C'^'scaia™"^ 
vs. m^, which extends also to light WIMP masses. The 
DAMA Collaboration analysis of Ref. [i^ takes into ac- 
count the same variability in galactic halo models of Ref. 
[i(tI |. which is also used here. It is not possible to make 
direct comparison among the DAMA allowed region and 
the upper limits we obtain here for CDMS and ZEPLIN, 
since the DAMA region is the convolution obtained af- 
ter varying all the galactic halo models, while the re- 
sults presented here refer to single halo model. A proper 
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comparison between different experimental results can be 
made only at a fixed galactic halo model. Notice that a 
convolution of our results would be just the upper curve 
(model Bl) of Figd 

As for the comparison between the light ncutralinos 
of non-universal gaugino models and the DAMA allowed 
region, we comment by reminding that these light ncu- 
tralinos are totally compatible with the allowed DAMA 
region, as we showed in Ref. 0: they could in fact ex- 
plain the annual modulation effect. 
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